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We have investigated the adsorption of Li on graphene oxide using density functional theory. We 
show a novel and simple approach to achieve a positive lithiation potential on epoxy and hydroxyl 
functionalized graphene, compared to the negative lithiation potential that has been found on pres- 
tine graphene. We included the van der Waals correction into the calculation so as to get better 
picture of weak interactions. A positive lithiation potential suggests a favorable adsorption of Li 
on graphene oxide sheets that can lead to an increase in the specihc capacity, which in turn can 
be used as an anode material in Li-batteries. We find a high specific capacity of ~860 mAhg“^by 
functionalizing the graphene sheet. This capacity is higher than the previously reported capacities 
that were achieved on graphene with high concentration of Stone-Wales (75%) and divacancy (16%) 
defects. Creating such high density of defects can make the entire system energetically unstable, 
whereas graphene oxide is a naturally occurring substance. 


I. INTRODUCTION 

Energy storage is a major challenge in present times to 
meet the demand of consistent power supply, both in the 
portable devices like laptops, cameras, cellphone, etc., 
and in the bigger devices like electric vehicles, backup 
inverters, etc. Additionally, increase in the cost of fuel 
and increase in green house gases have given rise to the 
hybrid electric vehicle technology with battery supple¬ 
ment, known as plugin hybrid electric vehicle (PHEV).^ 
For several years, PHEVs and many other appliances. 
Lithium batteries have been considered to be optimum 
candidates for rechargeable battery technology. 

Although Li metal has a high theoretical capacity of 
~3842 mAhg“^, rechargeable lithium batteries cannot 
use the bulk lithium metal due to safety issues, den¬ 
dritic growth, rapid decay and small recycling time.^^El 
To overcome these problems many other materials and 
structure have been investigated, e.g., metal oxide such 
as LiCo02, LiFeP 04 and LiMn02,^® where Li atoms 
can be stored in layers within these metal oxides. How¬ 
ever, these materials are found to have very small prac¬ 
tical capacities of 140 mAhg“^, 170 mAhg“^and 119 
mAhg“^, respectively. Whereas, graphite with theoret¬ 
ical capacity of ~372 mAhg“^is used as conventional 
anode. Recently, significant work is being done to en¬ 
hance the capacities of anodes as well as cathodes us¬ 
ing the carbon based materialslike a composite an¬ 
ode based on carbon and silicon with a specihc capacity 
of ~1000 mAh g~^has been proposed.® Similarly, a few 
more work^^E^ also shows a hike in the capacity of the Li 
storage on graphene porous networks with Stone-Wales 
(SW) and divacancy (DV) defects. They reported an 
experimental capacity of 850 mAhg“^, whereas theoreti¬ 
cal calculation predict capacities of 698 mAhg“^and 590 
mAhg“^with 75% SW and 16% DV defects, respectively. 
Although even higher capacities with 100% SW and 25% 
DV defects are predicted from the calculation, creating 
graphene sheets with such a high percentage of defects 


is difficult and energetically unfavourable. For example, 
formation energies of SW defect is ^5 eVEStH and it is 
formed by rapid quenching from high temperature under 
the radiation of 90-100 eV of electron energy. Similarly, 
formation energies of single- and di-vacancies in graphene 
sheet are ~7.5 eV and ~8 eV, respectively.E^tHl Although 
DV are thermodynamically more stable than single va¬ 
cancy at high temperature, these vacancies are practi¬ 
cally immobile .E^ Additionally, removal of a large num¬ 
ber of atoms can lead to the bending and wrapping of 
the sheet destroying its electronic properties. To over¬ 
come the problem, as well as taking the advantage of the 
fact that graphene sheet with non-stoichiometric nature 
can increase the lithiation potential (LP), we investigated 
the lithium adsorption over graphene oxide (GO). Inter¬ 
estingly, we find a positive LP on these functionalized 
graphene sheets suggesting a favourable adsorption of Li 
on GO leading to the specific capacity of ^860 mAhg“^. 
The most notable aspect of our proposed system is that 
it is simple and cost effective to prepare. Functionalized 
graphene can be created by simply treating it with strong 
acid /basd22H^ that results in functional groups (epoxy (- 
0-) and hydroxyl (-OH)) attached to the sheet. For GO, 
one does not have to go through the reduction process 
of removing oxygen that naturally gets adsorbed on it 
or one does not have to create any high energy defects. 
Functionalized graphene shows a great potential for Li 
storage for making high capacity batteries. A specihc 
capacity of ~860 mAhg“^ obtained by us is greater than 
that of 698 mAhg“^and 590 mAhg“^obtained with 75% 
SW and 16% DV of defects. Thus, GO is not only easy 
to prepare, it also gives very high specihc capacities and 
could emerge as an ideal material for Li-batteries. 


II. METHODS 

Density functional theory (DFT) calculations are per¬ 
formed in the pseudopotential formalism using Quantum 
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EspresscP^ code. Due to the presence of lithium metal, a 
higher plane-wave cutoffs of 60 Ry is used. Exchange in¬ 
teractions were treated within GGA with Perdew-Burke- 
Ernzerhof (PBE^^ functional form. We did include spin 
polarization in our calculations but found that there is 
negligible magnetic moment in the system. As van-der- 
waal’s interactions play an important role in binding 
of Lithium ions over the graphene sheet we hav e in- 

cluded dispersion correction in our DFT calculations.^^HMI 

We choose a computationally cost-effective “DFT-D2” 
method of Grimme,!^ for VdW interactions. For (1x1) 
unit cell of graphene, we use a Monkhorst-pack k-mesh 
grid of 18 X 18 X 1 and the corresponding k-mesh for big¬ 
ger supercells. All the configurations have been relaxed 
with a force tolerance of 10“^ Ry/bohr. To avoid inter¬ 
actions between images perpendicular to the plane of the 
sheet, a vacuum of 21 A is used. Using these parameters, 
we get C-G bond length of 1.42 A for pristine graphene 
that agrees well with the experiments.ISll 

We calculate the adsorption energy and the lithiation 
potential as a function of Li concentration in pristine and 
functionalized graphene. LP can be determined from the 
Nernst equation using following equation: 


LP = 


zF ’ 


( 1 ) 


where z is charge on the Li ion and AG/ = AF/ -|- 
TASf — PAVf is the change in Gibb’s free energy. At 
room temperature, the entropic effects (TASf) and vol¬ 
ume effects (PAVf) will be very small comparecP^to the 
formation energy AEf and can be neglected. AEf is de¬ 
fined as AE^ — \Ei^i^-\-suhs {pEj^i -t- where 

Ehir^+subs is the total energy of the lithiated GO, En is 
the total energy of a single lithium atom in body centred 
cubic (BGG) structure, and Esubs is the total energy of 
substrate GO. On the other side, the AE of Li is calcu¬ 
lated by using E^i for a gas phase Li atom instead of the 
metallic Li in the BGG structure. 


III. RESULTS AND DISCUSSION 
A. Adsorption of Li on graphene 

We have looked at the adsorption of Li ion on graphene 
for various Li concentrations by either varying the cell 
size or keeping the cell size fixed and varying the number 
of Li atoms in the cell. Figure shows the adsorption 
energies as a function of Li concentration from 0.01% to 
0.10% by varying the cell size from 3x3 to 6x6, while 
the inset shows the variation of Li concentration from 
0.11% to 0.33% by varying number of Li atoms in a fixed 
cell size 2x2. First of all, notice that adsorption en¬ 
ergy (AE) of Li increases with decrease in Li concentra¬ 
tion suggesting a repulsive interaction between Li atoms. 
AE at hollow site is the highest of all the other high 
symmetry sites (onsite, hollow and bridge), in agreement 



FIG. 1. Adsorption energies versus Li concentration from 
0.01% to 0.1% on graphene by varying cell size for hollow and 
onsite positions. The inset shows adsorption energies versus 
Li concentrations from 0.1% to 0.33% in a fixed cell size of 
2 X 2 at hollow site. 


with previous literature.^^^^^ We find that including van- 
der-waal’s (VdW) increases the adsorption of Li by 0.35 
eV irrespective of the Li concentration and the disper¬ 
sion contribution to the AE comes out to be 0.48 eV, 
in good agreement with 0.45 eV found by Allouche et 
alW^ Even though the AE changes by a fixed amount 
of ~0.35 eV for all the Li concentrations, this difference 
is expected to change when we functionalize graphene 
with epoxy/hydroxyl groups, as these functional groups 
will change the charge distribution in the system. Hence, 
just by taking the VdW contribution separately will not 
give you the exact information about the AE and the 
proper inclusion of “DFT-D2” inside the scf calculation 
is important. 


B. Lithiation of Graphene Oxide 

Here we investigate GO as an anode material for Li- 
batteries. It is well known that pristine graphene is 
not a suitable candidate for th e Li adsorption due to 
its negative lithiation potential,and hence it has 
a limitation to intercalate only between the spaces of 
graphene planes (with a theoretical capacity of ^372 
mAhg“^). Recently i t has b een shown that by introduc¬ 
ing defects/vacancies functionalizin^AMI the 

graphene sheet, a positive LP can be achieved, and hence 
one can enhance the battery capacity. 

The structure of GO is still a puzzle due to its non- 
stoichiometric nature. Most recent NMR experiments^ 
suggest that the Garbon:Epoxy:Hydroxyl content on 
graphene can be varied by changing the environmental 
conditions. Gommonly, the G:0 ratio in GO is found to 
be 2:1, however after reduction process, this ratio can 
be increased to as high a value as In our 

work, we use C:0 ratio as 18:1 and 9:1. Under 0-rich 
environment, the GO contains epoxy functional group (- 
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O- attached to bridge site), while in H-rich environment, 
-OH group is formed (attached to carbon atom with ori¬ 
entation along hollow site) as shown in Figure (I) and 
(II), respectively. We looked at the Lithiation poten¬ 
tial on three types of functionalized graphene: (I) epoxy 
covered, (II) -OH covered, and (HI) both, epoxy and - 
OH covered as shown in Figure [^. These calculation 
are performed in a 3x3 supercell. Numbers in the Fig¬ 
ure]^ show the positions of Li atoms with respect to the 
functional gronp. Arrows in Figure]^ (HI) indicate 
possible -OH positions and orientations with respect to 
the position of the epoxy group in same side as the epoxy 
group. 
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FIG. 2. Models of GO with three different functional groups: 
(I) Epoxy (II) Hydroxyl and (III) Epoxy-|-Hydroxyl, attached 
to graphene sheet. Squares with numbers show the adsorp¬ 
tion sites of Li. In model(III) blue arrows show the possible 
configurations (a-h) with different orientations on the same 
side of epoxy functional. 

For the case of single Li atom adsorption on the oppo¬ 
site side of epoxy group (site-3,4,5), we get the highest 
LP of 0.03 eV at site-3 having specific capacity of 115 
mAhg“^. For the adsorption at site-1 and 2 on the same 
side as epoxy group, we get a LP of 1.44 eV and 0.31 
eV, respectively. This high LP of 1.44 eV at site-1 is 
due to the binding of Li atom to the Oxygen atom caus¬ 
ing it to break the double bond with the graphene sheet 
(see Figure]^ (H)). A high LP of 1.66 eV was also found 
by Stournara et al.^^ for a similar configuration of Li 
bonded to the oxygen atom. This discrepancy of 0.22 
eV in LP could be due to the different cell sizes used in 
our calculations. At site-2, the Li atom remains at the 
hollow site and hence has a lower LP as expected. We 
also looked at adsorption of Li atom at site-3,4,5 on the 
same side as epoxy and found that in all three cases, Li 
atom gets bonded with the oxygen atom resulting in 1.44 
eV of LP(same as site-1 conf. Figure[^(H)), with 0-Li 
bond getting oriented along the hollow site. Thus, LP 
gets greatly affected whether the Li atom is deposited on 
the same side or the opposite side of the epoxy group in 


the sheet. For the case of two Li atoms in cell at site-l&2 
and site-l&4 placed along the side of epoxy, we obtain a 
LP of 1.11 eV and 2.01 eV, respectively. For site-l&2 one 
Li atom makes bond to oxygen and another Li atom re¬ 
maining at site-2, while for site-l&4 Li atoms come near 
to oxygen with the position of oxygen at onsite position 
shown in Figure]^ (HI). 



30000000 3 0 0 0 0 0 


(IV) (VI 



FIG. 3. Gonfigurations (I)-(IV) show different forms of oxygen 
obtained on GO sheet in presence of Li. Configuration (V) 
and (VI) show side and top view of the highest capacity (860 
mAhg~^) configuration with positive lithiation potential. 

In the case hydroxyl functionalized graphene sheet, we 
get the highest LP of 0.86 eV at site-6 and 7 with Li atom 
deposited on the opposite side of the -OH group. Thus, 
LP in the hydroxyl functionalized graphene is higher than 
that of the epoxy functionalized graphene. With two and 
three Li atoms adsorbed on the opposite side, LP of 0.32 
eV and 0.22 eV (specific capacity 230 mAhg“^and 345 
mAhg“^), respectively are obtained. In the extreme case 
of fully covered (nine Li atoms per -OH group) GO a 
LP of -0.20 eV is found. However when Li is placed on 
the same side as the -OH group, it gets strongly bonded 
to the -OH group and the whole molecule (Li-OH) gets 
detached from graphene sheet. 

We then investigated the LP when both groups (epoxy 
and -OH) are presents in the graphene sheet and placed 
on the same side. The possible orientations of -OH with 
respect to epoxy are shown in Figure[^(HI). Out of (a)- 
(h) configurations, the lowest energy configurations are 
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FIG. 4. Lithiation potential as a function of capacity for 
epoxy, hydroxyl, and both (-0 & -OH) functional group s at- 
tached to the graphene sheet and compared with others.^^^ 
Squares represent LP with functional groups, while circles 
represent data from previous literation with SW and DV de¬ 
fects in graphene sheet. A highest capacity of ~860 mAhg“^is 
obtained on GO with both epoxy and hydroxyl groups at¬ 
tached. 

found to be (e) and (h). Lithium atoms are place on the 
opposite side of the functional groups. LP for one Li in 
3x3 supercell (107 mAhg“^) is 0.76 eV, which is 0.1 eV 
less than that obtained with only -OH group. However, 
as we gradually increase the concentration of Li, LP in¬ 
creases compared to that with only OH group, as shown 
in Figure]^ (compare red squares with blue squares), 
e.g., LP at 315 mAhg“^is 0.19 eV greater than that with 
only OH group. The close analysis of the structure shows 
that, above 315 mAhg“^capacities, epoxy group breaks 
one of its bond and oxygen gets attached to a single C 
atom at onsite position in the sheet. We observe a simi¬ 
lar behaviour when epoxy contained graphene sheet was 
fully covered with Li on the opposite side (average LP = 
-0.19 eV) as shown in Figure]^ (IV). This new behaviour 
of oxygen atom replicates a situation similar to that of 
the hydroxyl group that shows a higher LP compared 
to the epoxy group. This breaking of the bond of oxy¬ 
gen with one of the carbon atoms could be caused by the 
stretching of the graphene sheet in presence of the metal¬ 
lic Li layer on the opposite side. Charge redistribution in 
presence of Li atoms (particularly the 2s orbital contribu¬ 
tion) and the singly bonded oxygen atom with graphene 
sheet at onsite position could be important factors for 
the boost in the LP with higher capacity. From Table [^, 
notice that the configuration HI(e) dominates the bind¬ 
ing upto the capacity of 645 mAhg“^with LP of 0.17 eV, 
while from 752 mAhg“^ capacity onwards configuration 
HI(h) gives the highest capacity ^860 mAhg“^with LP 
of 0.03 eV. Thus, we successfully achieve a positive LP 
upto the specific capacity of ^860 mAhg“^just by plac¬ 
ing the Lithium atoms only on one side of the sheet as 


shown in Figure]^ (V) & (VI). 

TABLE I. Lithiation potential with varying capacity for conf. 
Ill(e) and Ill(h) in 3x3 supercell. 


Gapacity 

Conf. Ill(e) 

Conf. Ill(h) # of Li 

(mAhg-i) 

eV 

eV 

atoms 

107 

6.76 

6.54 

1 

215 

0.29 

0.26 

2 

322 

0.41 

0.28 

3 

430 

0.24 

0.20 

4 

537 

0.26 

0.12 

5 

645 

0.17 

0.16 

6 

752 

0.08 

0.13 

7 

860 

0.02 

0.03 

8 

967 

-0.09 

-0.02 

9 


IV. CONCLUSIONS 

Lithium adsorption and lithiation potential on 
graphene oxide has been investigated using DFT. We find 
that epoxy and hydroxyl functionalized graphene sheet 
provides the highest specific capacity of ~860 mAhg“^. 
We have achieved higher capacities than previously re¬ 
ported capacitieJi^ of 698 mAhg“^and 590 mAhg“^with 
75% SW and 16% DV defects. Although, even higher ca¬ 
pacities of 1186 mAhg“^and 1029 mAhg“^were obtained 
by Dutta et a/.,E^with 100% SW and 18.75% (or higher) 
DV defects, respectively. We believe that creating such 
high density of defects would make the entire system en¬ 
ergetically unstable leading to bending and wrapping of 
the 2D sheet. In this paper, we have shown the hike in 
Li storage capacity by attaching the functional group on 
one side of the sheet and Li on the other side. How¬ 
ever, it would be possible to further increase the storage 
capacity by increasing the density of functional groups 
and Li atoms, Le., by playing around with number of Li 
atoms per epoxy and/or per hydroxyl groups. We have 
also demonstrated that London’s dispersion contribution 
plays an important role in binding of Li on graphene and 
GO sheets. Oxygen bonding and its position in GO plays 
an important role in raising the specific capacity of GO. 
Thus, based on the first principle studies, we propose an 
efficient way to increase the capacity on a novel low-cost 
system (GO). Our results provide an insight to find a 
better anode material for rechargeable Li batteries with¬ 
out damaging or drastically modifying the properties of 
graphene. 
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